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Strongly-interacting ultra-cold atoms in tight-binding optical lattice potentials provide an ideal 
platform to realize the fundamental Hubbard model. Here, after outlining the elementary single 
particle solution, we review and expand our recent work on complete characterization of the bound 
and scattering states of two and three bosonic atoms in a one-dimensional optical lattice. In the 
case of two atoms, there is a family of interaction-bound "dimer" states of co-localized particles that 
exists invariantly for either attractive or repulsive on-site interaction, with the energy below or above 
the two-particle scattering continuum, respectively. Adding then the third particle — "monomer" — 
we find that, apart from the simple strongly-bound "trimer" corresponding to all three particles 
occupying the same lattice site, there are two peculiar families of weakly-bound trimers with energies 
below and above the monomer-dimer scattering continuum, the corresponding binding mechanism 
being an effective particle exchange interaction. 

I. INTRODUCTION 

Among the tight-binding lattice models of condensed matter physics [J Q , the Hubbard modelQ plays a funda- 
mental role. It describes particle tunneling between adjacent lattice sites as well as short range (contact) interaction 
between the particles on the same lattice site. Despite apparent simplicity, this model is very rich in significance 
and implications for the many body physics on a lattice [J|. This is perhaps most profoundly manifested with nu- 
merous important experimental and theoretical achievements with cold neutral atoms trapped in deep optical lattice 
potentials wherein the Hubbard model is being realized with unprecedented accuracy. 

A remarkable Hubbard model phenomenon is the existence of stable repulsively-bound pairs of atoms in an optical 
lattice, as was experimentally demonstrated in Ref. [8| . This seminal achievement has led to several theoretical studies 
of the properties of interaction-bound atom pairs — "dimers" — in periodic potentials [9HL9|. 

In a one-dimensional (ID) tight-binding lattice, two bosons can form a bound dimer[8Hl3|] for any finite strength 
of the on-site interaction, be it an attraction or a repulsion. Next level of complexity corresponds to three bosons, 
which obviously can form a strongly-bound "trimer" with all three particles occupying the same lattice site. For 
large enough on-site interaction strength, however, there are two more kinds of weakly bound trimers 20]. The 
corresponding binding mechanism turns out to be an effective particle exchange interaction between the dimer and 
the third particle "monomer" leading to symmetric and antisymmetric trimer states with energies slightly above 
and below the continuum of scattering states of (asymptotically) free dimer and monomer. 

Below, after introducing the Bose-Hubbard model and outlining its elementary single particle solution, we review 
and expand our recent work [l2l. [20j| on the bound and scattering states of two and three bosonic atoms in a ID optical 
lattice. 

II. THE MODEL 

Cold bosonic particles in a ID tight-binding lattice can be accurately described by the (second-quantized) Hubbard 
Hamiltonian [3|-|7[ 

j 3 

where M (bj) is the particles creation (annihilation) operator and hj = btbj the number operator at jth lattice site; 
J(> 0) is the inter-site tunneling, or hopping, rate; and U is the on-site interaction, which can be attractive or 
repulsive. 

a. Single particle solution: Denoting by | j) the state with a single particle at the jth lattice site, the Hamiltonian 
reduces to 

H^ = -Jj2(\j)(j + M + li + i)01)- (2) 
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FIG. 1: Single particle Bloch band in a ID tight-binding lattice. 

Expanding the single-particle state vector as — J2j "00) \j) , the stationary Schrodinger equation H^ 1 ' \ip) = 
E^ \ip) leads to the difference equation 

- J + f ) + <Kj - 1)] = (3) 

which is satisfied by the (discrete) plane wave ansatz = V'gO') = cx P(^j) f° r the wave function with quasimo- 
mentum q £ restricted to the first Brillouin zone f2 = [— 7r,7r], while for the corresponding eigenenergy we obtain 

Eq = e(q) = — 2Jcos(q). The single-particle energy (Bloch) band, shown in Fig.Q] has, therefore, a width of 4 J. 

Having reviewed the trivial single-particle solution in the following Sections we present complete solutions of 
the much richer two- and three-body problems in a ID lattice. 




III. TWO PARTICLES IN A LATTICE 



Although the Bose-Hubbard Hamiltonian ([T]) corresponds to identical bosons, which is our main concern here, it is 
nevertheless instructive to tackle the more general problem of two distinguishable particles in a lattice (Tlj . a simple 
limit of which produces the solution for two indistinguishable bosons. 

We thus consider two particles A and B having, in general, different hopping rates J a and Jb, respectively. The 
Hamiltonian describing their dynamics is the non-symmetrized (first-quantized) version of the Hubbard Hamiltonian 
(TTJ) acting on the two particle subspace: 

ff (2) = -J A ^(|jU>0U + l| + \ja + 1)(3a\)-JbJ2(\3b)(3b + 1\ + \jb + 1)Ub\) + U £ \ Ja ,Jb)(3a,3b\ , (4) 

J A JB j A —JB 

where ja and js denote the positions of particles A and B, respectively. The eigenstates of can be expanded as 
Y$>) = J2j A j B ^Ua, Jb) \j A ,jB), so that the Schrodinger equation \^>) = E^ |^) leads to the difference equation 

- Ja^Ua + 1,3b) + *(M - 1,3b)] ~ Jb^UaJb + 1) + *(j A ,3B - 1)] + U5 jA , jB ^(j A ,j B ) = E^(j A ,j B )- (5) 

In order to solve the problem analytically, we need to transform the two-body difference equation ([5]) into a "one- 
body" problem. To that end, we define the center of mass jr = h (j A + 3 b) and relative j r = 3a ~ 3b coordinates and 
use for the two-particle wave function the separation ansatz [l5ll20| 

*Va,3b) = e^e-Wri'lucUr), (6) 

where 

tan(fe) = i A 7 J f tan(if/2), (7) 
J a + Jb 

with K £ Q the center-of-mas s q uasimomentum. Note that for J a = Jb~J the separation ansatz (jSJ) reduces to 
that for identical particles [l2l [l4j . The resulting recursion relation now reads 

- J K [i>K(jr + 1) + lP K (jr - 1)] + US 3r ^J K ( Jr ) = Efi, K {j r ), (8) 
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where Jk = \J J\ + J\ + 2JaJb cos(K) is the collective hopping rate [HHHI, which in the case of identical particles, 
J A = Jb = J, reduces to the standard [1, H, [H, [H| expression Jk = 2Jcos(if/2). 

Equation flHJ admits two kinds of solutions, corresponding to the scattering states of asymptotically free particles 
and to the two-particle bound, or dimer, states. 

A. Scattering states 

Since the interaction between the particles is governed by a short-range — in this case a contact, U5j rj o — potential, 
its action amounts to a unitary phase shift (see below), while the spectrum of such solutions is given by the sum of 
the spectra for two free particles A and B with momenta qA = K/2 + k and gs = K/2 — k: 

E K, k = 4a + 4b = - 2J K COS ( fc )> ( 9 ) 

which spans the interval E K ' k G [— 2Jk, 2Jk\- The corresponding symmetric scattering wave functions, for sin(fe) ^ 0, 
are given by 

1pK,k{jr) = C0s(fc|j r | + S K ,k), (10) 

which, upon substitution into Eq. ([8]). yields for the scattering phase shift 5k,)c, 

tfw(5 K ,k) = - T -r-TTs - (11) 

2 Jk sin(re) 

For U — > we have non- interacting particles ipK,k{jr) — cos(fc|j r |), while for U/ Jk ±oo we obtain the fermionized 
solution ipK,k(jr) = sm (fc|jr|), whereby the two particles never occupy the same lattice site, t^K,k(0) = 0. 

When sin(fc) = 0, the scattering wave functions (flU)) are no longer valid, and the lattice generalization of the 
continuum zero-energy solution of the Schrodinger equation apply [Tol ll5|| . At the bottom and the top of the scattering 
band, E^' T = ^2J K [k = 0,tt), the corresponding solutions ipK,^ have the form 

1> K ,-(jr) = 1-^L, (12a) 
a-K- 

j>K,+ Ur) - (-1) J " U - , (12b) 



where clk,^, with o>k,- — —ax.+, are the scattering lengths, which are calculated by substituting Eqs. (|12p into 
Eq. ([8]), resulting in <ik.- = —2Jk/U (in units of the lattice constant). Note that, since the scattering lengths are 
finite for any U ^ 0, there exists no "zero-energy" resonance in this model. Therefore, for every value of the total 
quasimomentum K there is only one bound state, as we shall see below. 

To illustrate the foregoing discussion, in Fig. [5] we show the energy spectrum © for two identical bosons, J a = 
Jb = J, and the corresponding density of states defined via 

L dk _ L 1 
P{E ' K) ~2^dE ~2^^\U^(K/2)Y-E^ (13) 

with L a quantization length. As seen, the density of scattering states, i.e., the number of states per unit interval of 
energy, is lowest in the middle of the band, E ~ (and K ~ 0), while for a given quasimomentum K of the center of 

( 2) 

mass motion of the two particles, p(E, K) increases rapidly as energy E = E K k approaches its maximal and minimal 
values ±2Jk- 

B. Bound states 

The two-body bound states of Hamiltonian Q are the solutions of Eq. (|SJ) yielding normalizable relative coordinate 
wave function, J^j I^K(ir)| 2 < oo, with the corresponding energy below (for U < 0) or above (for U > 0) the 
scattering continuum at each value of the total quasimomentum K. Introducing into Eq. ((8} the exponential ansatz 
ipK(jr) oc a K r)[ yields the bound state energy E^ = —Jk{olk + 1/or-), with given by 
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FIG. 2: Energies E K versus the center-of-mass quasimomentum K for a pair of identical bosons in a ID lattice. The continuum 
spectrum corresponds to energies Q of the scattering states, with the shading proportional to the density of states (|13[) . The 
(blue) line below and the (red) line above the scattering band are, respectively, the energies (|14|l of the attractively-bound 
dimer with U = — 5J and the repulsively-bound dimer with U = 5J. 

We see that olk is real and also \<Xk\ < 1 for all U ^ 0. Hence, the relative coordinate wave function ipK(jr) is 
normalizable, and the energy satisfies \E K \ > 2Jk for all K. Explicitly, the energy and the normalized wave function 
for the bound dimer are given by 

E%> = sgn(U)^U 2 + U 2 K , (14) 

<M» = ^=2= (U K - sgn(U)Ju 2 K + l) brl . (15) 
V "if + 1 

In Fig. [2] we show the energies E K for a pair of identical bosons, Ja = Jb= J, interacting via on-site attractive, 
U < 0, or repulsive, U > 0, potential. The binding energies E 1 ^ are defined with respect to the edges of the scattering 
band ([9]). Note that in the case of repulsive interaction, U > (o.k < 0), the sign of the wave function ((T5|) alternates 
between the neighboring sites j r . Remarkably, when \K\ = tt, and thereby Jk = 0, the relative coordinate wave 
function tpK(jr) is completely localized at j r = for any [7^0. 

Clearly, for a given value of the dimer quasimomentum K, the stronger is the on-site interaction \U\, the smaller is 
the extent of the wave function ipi<(jr), meaning that the constituent particles are stronger co-localized. For \U\ 3> J, 
the dimer energy in Eq. (|14[) can be approximated as 

Ef ~2J^cos(K) 7 (16) 

where the first term E ^ = [U — 2J^] represents the dimer "internal energy", while the second term e^ 2 '{K) = 
—2,j( 2 ) cos(K) is the kinetic energy of a dimer with quasimomentum K and an effective tunnelling rate J*- 2 * 1 = —2J 2 /U. 

IV. THREE PARTICLES IN A LATTICE 

Building on the solution of the two-body problem, in this section we consider three bosonic atoms is an optical 
lattice described by Hamiltonian (P). 

A. Spectrum of scattering states 

With the expertise gained from the previous section, we can readily deduce that in the case of three particles there 
are two distinct scattering continua (see Fig. [3] left panel). The first is the three-body scattering continuum of three 
(asymptotically) free particles, with the energy given by the sum of single-particle bands, E c3 = e(fei) + e(k 2 ) + e(K — 
k\ — k 2 ) = e(ki, &2, K — k\ — k%), where each particle quasimomentum kj € f2 is in the first Brillouin zone = [—tt, tt] 
and K = k\ + k% + ^3 (mod 2tt) is the total quasimomentum. The second is the two-body scattering continuum of a 
bound pair (dimer) and a free particle (monomer), with energy E c2 — sgn(U)y/U 2 + [4 J cos(Q/2)] 2 — 2Jcos(K — Q), 
where the first term is the energy of a dimer, Eq. (|14p. with quasimomentum Q. 
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B. Bound states 



We seek the bound states j*) of three bosons in momentum representation, 



I*) 



(2tt) 3 / 2 



dkidk 2 dk 3 k 2 , k 3 ) \ki, k 2 , £3), 



(17) 



where the wave function /c 2 , /C3) is symmetric with respect to exchange of any pair of particles. From the 

stationary Schrodinger equation H |^) = _£7 1 , using the conservation of total quasimomentum K, we obtain [2l| 



*(fc 1 ,fc 2 ,fc 3 ) = 



M(fe x ) + M(k 2 ) + M(k 3 ) 



e(fci,fc 2 ,fc 3 ) 

where functions M(k) satisfy the ID Mattis integral equation plj 



E 



M(k)[l + I E (k)] 



with lE(k) being a generalized Watson integral [22 



U 



dq 



M(q) 



e(k,q,K - k - q) - E' 



(18) 



(19) 



I E (k) 



U 
2^ 



dq 



sga[E-e(k)]U 



e(k,q,K-k-q)-E ^[E - e(k)} 2 - 16 J 2 cos 2 [(K - k)/2] 



Equation (fH?)) can be cast as a homogeneous Fredholm equation of the second kind with eigenvalue A = 1. Hence, for 
a given U/J and fixed K, it is a nonlinear equation for energy E, which we solve numerically. 

The full three-body spectrum of Hamiltonian ([1} is shown in Fig. |3]left panel, with the bound state energies denoted 
by Ef,, E a \ and E a2 . For concreteness, here we consider attractive interaction, U < 0, but note that our results equally 
apply to the case of repulsive interaction, U > [HI, [2(| . 

In complete analogy with the two-body problem, the Bosc- Hubbard Hamiltonian ((T|) with \U\/J 1 has a very 
narrow band of on-site bound states, corresponding to three tightly bound bosons co-localized on the same lattice 
site, with energies Eb « 3U far from both scattering continua 21] . But as also seen in Fig. [31 in a ID lattice, bosons 
can form two new kinds of three-body bound states whose energies E a \ and E a2 lie below and above the two-body 
continuum E c2 . Some properties of these states can be deduced by energy considerations. First, these are not on-site 
bound states, since their energies -E a i(2) — U + 0(J) are far from 3U. Next, their binding energies, with respect to the 
E c2 band, are Ebx{2) T>^/2, which suggests that these are off-site weakly-bound states of a dimer and a monomer. 
Note that the state above the two-body continuum is bound stronger than the state below the continuum. Finally, 
they are not Efimov states which can exist only in 3D systems near two-boson resonances [2ll . |23j . 

As can be seen from Fig. [3]right panel, where we plot the binding energies Ebi(2) at quasimomentum K — 0, there 
are thresholds for the existence of full bands (K € [— 7r, tt]) of the off-site bound states. For the trimer below the 
two-body continuum, the binding energy vanishes when \U\ w 4J: at this critical value of U the trimer energy E a \ 
approaches the edge of the dimer-monomer scattering continuum E c2 = ~^/U 2 + 16 J 2 — 2 J. On the other hand, 
the K = trimer above the two-body continuum ceases to exist already for \U\ ~ 8.5J, since then its energy E a2 
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FIG. 3: Left: Full three-particle energy spectrum of Hamiltonian |T]) with U = — 10J, versus the total quasimomentum K. All 
bound states are obtained via exact numerical solution of Eq. (|19[) . Right: Binding energies Eb for the off-site (weakly-bound) 
trimers at K — versus the interaction strength U < 0. 
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approaches the bottom of the three-body continuum E c % = 3e(0) = — 6 J (the two continua, E c2 and E c ^, overlap for 
\U\ < 8J). Thus, at K — 0, the trimer state with energy E a2 starts to appear well in the strong interaction regime, 
while for larger K the threshold is smaller: \U\ ~ 4J for \K\ — > tt. 



C. Effective model 



Since for strong on-site interaction \U\/J 1 the dimer is essentially unbreakable, the off-site trimers with energies 
close to U must be bound due to a mechanism different from the on-site interaction alone. To identify such a 
mechanism, we derive an effective perturbative model, valid for \U\/ J > 8, describing two distinguishable, hard-core 
particles — the dimer and the monomer. To second order in the tunneling rate J, the effective Hamiltonian reads 

H eS = H x + H 2 + fltat, (20) 
ffi = -J$Z($t$ i+ i+H.c.) 



where 



describes the single monomer; 



H 2 = £^ ]T mj - J (2) $>tc i+1 + H.c. 

j 



is the Hamiltonian for a dimer (c/. Eq. (|16[) ). with at being the dimer creation (annihilation) operator and 



rh j = c -Cj the number operator at site j; and finally 



Hurt = mjh j±1 - wY,{c] +1 c 3 b]b J+1 + H.c. 
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describes effective interactions between the dimer and the monomer, including a weak nearest-neighbor interaction 
= —7J 2 /2U, and an exchange interaction with the rate W = 2 J. As we will see below, it is the exchange term 
that is responsible for the formation of the off-site trimers. 

In Fig. |4] right panel, we plot the spectrum of the effective Hamiltonian (f20|) . which contains two bound states with 
energies E a ± and E a2 below and above the two-body scattering continuum E c2 = £^ + e^'{Q) + e(K — Q). These 
dimer-monomer bound states are obtained using the Schrodinger equation for the two-body wave function ^(Q, k) 
in momentum space, which leads to the integral equation 

k) = -tojj* £ m + eW{ q ) + e{K- q )-E * (q > ^ (21) 

where K = Q + k is the total quasimomentum, V cos (Q,q) = [2V^ cos(q) — 4J cos(K — q)} cos(Q) and V s i n (Q,q) = 
V cos (Q,q) with cos — ¥ sin, while Ui 2 (— ¥ oo) is an artificial dimer-monomer on-site interaction imposing the hard- 
core condition on Hamiltonian (|20[) . Equation (|21j) reduces to a non-linear equation for the energy E solving which 
(numerically) we obtain E a i and E a2 . Comparison with the exact spectrum on the left panel of Fig. |4] reveals good 
agreement: the continuum spectra are indistinguishable, while the small but noticeable differences in the bound-state 
energies are associated with the internal structure of the dimer, not accounted for by the effective model, and they 
gradually disappear with increasing the on-site interaction strength U. 




FIG. 4: Left: Magnified part of the exact three-particle spectrum of Fig. [3] Right: Dimer-monomer spectrum of the effective 
Hamiltonian (|20p with the two bound states obtained via numerical solution of Eq. (|21|l . 
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1. Analytic solutions 



There are two important cases, corresponding to the maximal (K — n) and minimal (K = 0) total quasimo- 
mentum K, for which the bound and scattering states of the effective Hamiltonian (|2"U1) can be calculated ana- 
lytically employing the method of Sec. IIII1 To that end, we expand the two-particle eigenstates in coordinate 
basis 1^) = ~Y^j 1 J.j 2 ^{ji,ji)\jitj2l with the wave function in the form ^(ji,^) = e lK; > R e~ l P K: > r 'ipK{jr), where 
j R = -(ji + j 2 ), j r = ji — J2, with ji and j 2 being the lattice positions of the monomer and dimer, and 
tan(/?R-) = tan (K/2)[J — J^]/[J + J™]. For the relative coordinate wave function ^px(jr), imposing the hard- 
core condition ipK (0) = 0, we then obtain the difference relations 

Jk^k(±2) + W k ^k(tI) + [E-V {2) }^k(±1) = 0, 

JK[lpK(jr + l) + lpK(jr-l)]+Eip K (j r ) = 0, (22) 



with \j r \ > 1, J K = y/J 2 + J( 2 ) 2 + 2JJ( 2 ) cos(K), W K ee Wcos(K), and E = E - £ (2 \ 

a. Scattering solutions: Using the standard ansatz ip^jjjr) = cos (fe|j' r | + ^kI) ano - ^Kktir) = 
sgn(j r ) cos (fc|jV| + &Kk) w ith k the relative quasimomentum, for the corresponding phase shifts S^f, of the sym- 
metric (+) and antisymmetric (— ) scattering wave functions we obtain 



? (±) ^ _ Jk cos(2fc) + [E± Wk — V^] cos(fc) 
J K sin(2fc) + [E±W K - V( 2 )] sin(fc) 



. / —I — \ , tj K IjUO ZjM.' I J_/ _L_ V V K V >^WO \ rv i , 



with E — E C 2 — £^ = —2Jk cos(fc). Note that in the limit of \U\ — > oo, as the the nearest neighbour interaction 
and the dimer hopping tend to zero (Jk — > J), Eq. (f23)) holds for all K = k. The full scattering wave function is 
given by a superposition tpK,k(Jr) = ^V'i-fc(jr) + B ipj^iUr), which, upon expressing through incident, reflected and 
transmitted waves, 



, fe ife > +r K e-^ (j r < 0) 

^'H tjr e^ ( ir >0) 



>0)' 

leads to v4/£? = — e~^ SK > k ~ SK > h \ For the reflection r^-(fc) and transmission tx(k) amplitudes we then obtain tk^k — 
l\^ e 2lS K,k -j- e 2l5 A-,fcj ; arLC [ the transmission and reflection probabilities are given by Tj( = |£/f| 2 = sin 2 {S^\ — 8^\) 
and R K = \r K \ 2 = cos 2 

In Fig. [5] we plot Tx{k) for U/J = —10 at total quasimomenta if = and \K\ = 7r. The transmission spectra for 
the intermediate values of K lies in between the curves for K = and K — tt. We observe the maximum transmission 
in the vicinity of k = ±7r/2, where Tk(it/2) ranges from 50% to 80%. With increasing the interaction strength U, we 
find that the maximum transmission saturates at around 64% for all values of K, which should be contrasted with 
the results of Ref. [lj|. 

b. Bound solutions: The dimer-monomer bound states are obtained from Eq. (|22p using the exponential ansatz 
ipxUr > 0) oc o?k~ X and ip K (-jr) = ±?pK{jr), which yields 

a (±) = (24) 

for the symmetric (+) and antisymmetric (— ) wave function of the bound state (\c(k\ < 1), with the corresponding 
energy E al{2) = -J K [1 + (a£>) 2 ]/a£°. 

It is now easy to see that without the exchange interaction there would be no dimer-monomer bound states (for 
any K). Indeed, this hypothetic (W — 0) problem is exactly solvable for all K, and for two hard-core bosons with 



nearest-neighbour interaction there could be only one bound state 14 [ when \olk\ = \ Jk/V^ \ < 1, which cannot 
be satisfied in the range of validity (\U\/J > 8) of the effective model, Eq. (f2"0")l . The effective nearest-neighbour 
interaction is, however, responsible for the asymmetry in the binding energies -E_bi(2) = ^ai(2) T 2Jk of the exchange- 
bound trimers below and above the continuum E c2 (see Fig. [3] right panel). With increasing the on-site interaction 
U , the binding energies attain the universal limits E B1 ^ 2 ) ~> which remain valid for all K, since — > as 

\U/J\ — > oo. We thus conclude that the effective particle exchange interaction binds the dimer and the monomer into 
the weakly bound trimer. 
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FIG. 5: Transmission probability Tr-(/c) of a single particle with relative quasimomentum k through a bound dimer, for 
U/J = —10 and total quasimomenta K = (red solid line) and \K\ = n (blue dashed line). 

V. CONCLUSIONS 

In this paper, we have presented complete solutions for the one-, two- and three-body problems in a ID tight-binding 
lattice described by the Bose-Hubbaxd model. 

For the case of two bosons, we have derived the scattering states and the bound states of co-localized particles, 
which we termed as dimers. The corresponding binding mechanism is provided by the on-site interaction between 
the particles. Perhaps counterintuitively, these bound dimer states exist not only for attractive, but also for repulsive 
interactions Q. In free space, the repulsive interaction would inevitably lead to the pair separation, or dissociation, 
whereby the potential energy of the repulsion is converted into the kinetic energy of the free particles. In the lattice, 
however, the kinetic energy of each particle cannot take on arbitrary values, but is restricted to the values in the 
allowed Bloch band, which itself is bounded both from below and from above. Therefore, in the absence of energy 
dissipation, a pair of co-localized particles interacting even via repulsive potential is destined to stay bound together 
as a dimer, simply because there are no available free-particle energy states to which the dimer can dissociate. 

For three bosons, we have found three families of trimers — bound states of three particles. The first strongly-bound 
trimer, being an analog of the dimer, corresponds to all three particles occupying the same lattice site and bound 
by the on-site interaction. The other two families of trimers are weakly-bound with energies just below and above 
the two-body scattering continuum of a single particle (monomer) and an interaction-bound dimer. Intuitively, these 
trimer states correspond to symmetric and antisymmetric states of a dimer and monomer at the neighboring lattice 
sites interacting with each other via an effective (particle) exchange interaction, which is responsible for their binding. 

The phenomena discussed above are pertinent to the experiments with cold bosonic atoms in optical lattices [7|. 
Studying larger number of bosons in a lattice might reveal other exotic bound states, while longer range interactions, 
such as those between dipolar atoms [24j . or molecules p5| . will certainly play an important role in the formation of 
few-body bound states [lid HH • 
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